c      mainprogram mcube.f
c
c
       parameter (NX=21,NY=21,NZ=21)
       parameter (NGX=20,NGY=20,NGZ=20)
       parameter (D=10)
c
       dimension gc(3)
       dimension grc(3)
       dimension bc(3)
       dimension rc(3)
       dimension v(NX,NX,NX)
c
       gc(1)=1
       gc(2)=1
       gc(3)=1
c
       rc(1)=1
       rc(2)=0
       rc(3)=0
c
       grc(1)=0
       grc(2)=1
       grc(3)=0
c
       bc(1)=0
       bc(2)=0
       bc(3)=1
c
       call initvrml
c
       call viewpoint(NGX,NGY,NGZ,D)
       call dr_scale(gc,NGX,NGY,NGZ,D)
c
       call mk_data(v,NX,NY,NZ,D)
c
       tr=0.6
       th=NGZ*D
cc
       call mk_equ_face((th*1)/3,v,rc,tr,NX,NY,NZ,D)
       call mk_equ_face((th*2)/3,v,grc,tr,NX,NY,NZ,D)
       call mk_equ_face(th,v,bc,tr,NX,NY,NZ,D)
c
       stop
       end